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ABSTRACT 

The color of galaxies is a fundamental property, easily measured, that constrains models of galaxies and 
their evolution. Dust attenuation and star formation history (SFH) are the dominant factors affecting the color 
of galaxies. Here we explore the empirical relation between SFH, attenuation, and color for a wide range of 
galaxies, including early types. These galaxies have been observed by GALEX, SDSS, and Spitzer, allowing 
the construction of measures of dust attenuation from the ratio of infrared (IR) to ultraviolet(UV) flux and 
measures of SFH from the strength of the 4000A break. The empirical relation between these three quantities is 
compared to models that separately predict the effects of dust and SFH on color. This comparison demonstrates 
the quantitative consistency of these simple models with the data and hints at the power of multiwavelength 
data for constraining these models. The UV color is a strong constraint; we find that a Milky Way extinction 
curve is disfavored, and that the UV emission of galaxies with large 4000A break strengths is likely to arise 
from evolved populations. We perform fits to the relation between SFH, attenuation, and color. This relation 
links the production of starlight and its absorption by dust to the subsequent reemission of the absorbed light 
in the IR. Galaxy models that self-consistently treat dust absorption and emission as well as stellar populations 
will need to reproduce these fitted relations in the low-redshift universe. 

Subject headings: galaxies:fundamental parameters — galaxies:evolution — dust:extinction — ultravio- 
let:galaxies — infrared:galaxies 



1. INTRODUCTION 

Since the work of iTinsleyl (1 19681) the location of galaxies 
in color-color and color-magnitude diagrams has proven to 
be an important indicator of the stellar content of galaxies, 
just as the location of stars in these same diagrams indicates 
their stellar spectral type. This stellar content is in turn used 
to constrain the star formation histories (SFHs) of galaxies, 
which is an important step toward understanding the diversity 
of galaxies observed in the universe today. Color-color dia- 
grams are still used at high redshift (and low), where in ad- 
dition to the stell ar populations and dust attenuation, the red- 
shift is estimated ( Schweizer & Seitzer 1992; Mob asher et al.l 
l2004tlWuyts et al.ll2007l) . 

To obtain strong constraints on the constituent stellar popu- 
lations with color-color diagrams, it is important to know the 
properties of the dust attenuation in galaxies. Selective atten- 
uation can cause a large amount of reddening that can mask 
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the differences in color due to differences in stellar popula- 
tions. Such attenuation is oft en treated as a reddening 'vec- 
tor' in color-color space (e.g.. iKauffmann et al.ll2003bl) . The 
direction of the vector is given by the shape of the attenua- 
tion curve (often assumed to be that of lCalzetti et al.l (119941) ). 
while the length of the vector is inferred indirectly. Without 
independent knowledge of the length of the reddening vec- 
tor, the power of color-color diagrams to probe the constituent 
stellar populations of galaxies is reduced. 

The proper treatment of dust attenuation is of paramount 
importance when considering the constraints imposed by the 
color-magnit ude diagram (CMP) on theories of g alaxy for- 
mation (e.g., iNoeske et alj|2007t lLabbe et all [2007b . This is 
especially true since th ere is a strong correlatio n of atten- 
uation with luminos ity (Wang & Heckman 1996) and stellar 
mass (e.g., ISalim et al.l l2006)" These correlations can induce 
or alter trends between the derived SFH and stellar mass of 
galaxies. Similarly, dust attenuation also affects comparisons 
between the predictions of semi-analytic modelling and ob- 
servations of galaxies at any wavelength. 

With the Galaxy Evolution Explorer (GALEX) and Spitzer 
Space Telescope it is now possible to obtain ultraviolet (UV) 
and infrared (IR) fluxes for large samples of galaxies. The 
combination of the UV and IR flux constrains the amount of 
dust attenuation in galaxies (i.e. the length of the reddening 
vector). The UV is due to emission from young stars that is 
not attenuated by dust, while the IR measures the emission 
from young stars that has been absorbed by dust and reradi- 
ated. With this independent knowledge of the dust attenua- 
tion, which significantly affects the observed color of galax- 
ies, these colors can be used to more sensitively probe the 
stellar content of galaxies and test models of galaxy forma- 
tion. New models of galaxy spectra are being developed and 
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improved that treat dust attenuation, dust emission, and stel- 
lar por>ujatiojisin a cons istent manner, building from the work 
oflCharlot & Falll d2000l). These UV through IR models (e.g., 
ISilva et al.ll 1 998t Ijonssonll2006h promise to improve the treat- 
ment of attenuation and sharpen the constraints imposed by 
UV and IR observations on theories and semi-analytic mod- 
els o f galaxy formation. 

In I Johnson et al.l J2006) we presented the relation between 
broadband UV through near-IR colors, dust attenuation (as 
measured by the infrared excess, IRX), and SFH (as measured 
by the 4000A break) for a sample of galaxies observed with 
GALEX, the Sloan Digital Sky Survey (SDSS), and Spitzer. 
Here we present this same relation, with deeper UV data, for 
a larger sample of galaxies and a greater variety of broadband 
colors. We show that this relation serves as a useful diagnos- 
tic of galaxy properties, and an empirical constraint on mod- 
els of dust attenuation and star formation (SF). We compare 
our r esults to models of gala xy spectra dBruzual & Charlotl 
120031) and dust attenuation dWitt & Gordonl 120001) . In ad- 
dition, we consider where dwarf galaxies and IR- and UV- 
luminous galaxies fall in this diagram. This comparison of 
special classes of galaxies to a more representative sample 
places them in context, and aids the interpretation of the rela- 
tion. 

We present parameterized fits to this relation between color, 
the 4000A break strength, and IRX. These fits may be used to 
determine the dust attenuation in galaxies where only opti- 
cal spectra are available, even in the absence of measurable 
Ha and H/3 lines. In addition, these relations provide an em- 
pirical guide for self-consistent models of the UV through IR 
spectra of galaxies. 

2. DATA 

Our primary sample i s selected from the SDSS main galaxy 
spectroscopic sample (Straus set al.l l2002h . These galaxies 
have been imaged in five optical bands (ugriz) and observed 
spectroscopically from 3800A to 9000A (observed frame) in 
a 3"diameter aperture as part of the SDSS. 

In this paper we consider a subset of the SDSS spectro- 
scopic galaxy sample that has been observed by both GALEX 
and Spitzer. These galaxies are located in two different con- 
tiguous regions, the Lockman Hole and the Spitzer extra- 
galactic First Look Survey (FLS). Details of the observa- 
tions for each of these fields are given in Table Q] The 
Spitzer ob servations include four-band Infrared Array Cam- 
era (IRAC;SIiioItIH2004) imaging (3.6/im, 4.5/xm, 5.8/im, 
and 8/mi) as well as 24 an d 70 /im Multiban d Imaging Pho- 
tometer for Spitzer (MIPS; iRieke et all 2004) imaging. The 
GALEX observations include imaging at 1517A (far-UV, /) 
and 2261 A (near-UV, n). All quoted magnitudes are on the 
AB system and have been corrected for foreground extinction 
(which is especially low in the case of the Lockman Hole) 
according to the maps of lSchlegel et al.l (119981) . 

2.1. Optical Data and Derived Parameters: SDSS 

A description of the SP SS optical reductions is given in 
IStoughton etail (120021) and lAbazajian et all ((2004). Through- 
out this study we use the Petrosian magnitudes supplied by 
the SDSS Data Release 4 (DR4), as given in the photometric 
catalogs provided by the MPA/JHU group for SDSS studies 1 1 . 
Additional measurements of emission-line strengths and spec- 
troscopic index values have been made by the MPA/JHU 

1 1 http://www.mpa-garching.mpg.de/SDSS/ 



TABLE 1 
Observations 



Field Name 


Size 












(deg 2 ) 




f,n< 25 


m 2 4 < 19.5 




Lockman Hole 


~9 


872 


792 


819 


721 


FLS 


~3 


186 


147 


158 


118 



NOTE. — The median redshift in both fields is z = 0. 1 
a The number of SDSS/MPA galaxies that are observed at all / through 
24/^mwavelengths ( j|2.2t . b The number of observed SDSS/MPA galax- 
ies with fluxes brighter than the noted AB magnitude. The number of 
observed SDS S/MPA galaxies that pass the flux, redshift, and size cuts 
given in 92.41 

group using special-purpose code to improve continuum sub- 
traction, and are available from the MPA/JHU Web site. We 
have removed duplicate galaxies in these catalogs. 

Additional parameters including star-formation rates 
(SFRs) and stellar masses have also been derived and 
made publicly available for many of these galaxi es by 
iKauffmann et all d2003bl) and lBrinchmann et al.l (120041). The 
stellar mass estimates are obtained by iKauffman n et al.1 
d2003bl) from fits of a suite of stellar populution synthe- 
sis models to observed spectral indices (to obtain the dust 
free mass-to-light ratio) and colors (to estimate the atten- 
uation) . These stellar masses are consistent with stellar 
masses der ived using a simple magnitude and color prescrip- 
tion as in iBell & de J ong (2001), that are in turn consistent 
with limits from dynamical masses. The gas-phase metal- 
licities of many of the s ample galaxies have been deter- 
mined by Tre monti et al.l (120041) from emission-line fluxes. 
These metallicities are only available for galaxies with sig- 
nificant emission lines that are not classified as AGN (see 
IKauffmann et alJfeOOSal for the AGN classification criteria), 
and are only defined within the SDSS aperture. 

We only consider galaxies for whic h a mass has been de- 
termined by IKauffmann et all d2003bl) (here after SDSS/MPA 
galaxie s). Readers are encouraged to see Kauffman n et al.l 
(2003b) for a description of the optical selection. In total we 
consider 872 SDSS/MPA galaxies in the Lockman Hole that 
have been observed from the Far-UV through 24 /im. Of these 
galaxies, 866 have also been observed at 70 /im. In the FLS 
field we consider an additional 186 galaxies ob serve d in the 
UV (with sufficient UV effective exposure time, j32.21 i through 
24/im. 

2.2. IR and UV Data 

The Lockman Hole field has been observed by GALEX in 
26 0.°6 radius circular tiles. The effective exposure times of 
these tiles are ~ 10- 15ksin/ and ~ 10-60ks in n. In the FLS 
field the effective exposure times are highly variable from tile 
to tile, and we require exposure times > 5ks in both / and n. 
The UV photometry is taken from the pipeli ne-produced cat- 
alogs provided by the GALEX science team (Morrissey et al. 
120071) . For each SDSS/MPA galaxy we search the catalogs 
for objects within 3"of the SDSS position, choosing the clos- 
est object in the rare case of multiple matches. Since the 
GALEX fields overlap, some objects appear in two separate 
fields. We choose the tile with the longer effective exposure 
time (that takes into account sensitivity variations in the de- 
tector). We also exclude objects > 0.55 deg from the GALEX 
field center as the incidence of artifacts is much higher here. 
The UV photometry, the principal UV data used in this study, 
is computed in elliptical Kron apertures using a version of 
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SExtractor dBertin & Arnoutsl 1 1996t) modified to operate on 
low-background images. 

The IR data are from publicly available Spitzer imaging. In 
the Lockman Hole we use the reduced images provided by 
the Spitzer Wide-area Infrared Extragalactic Survey (SWIRE; 
lLonsdale et al.l 12003]) team 12 . In the FLS field we use the 
post-basic calibrated data (post-BCD) data provided by the 
Spitzer Science Center 13 . We have checked the consistency 
of these data sets by comparing the post-BCD Lockman Hole 
mosaics to the SWIRE-processed data in this field; the pho- 
tometry is consistent to better than 0.1 mag at the flux lev- 
els we consider, except at 70 /im (which we do not consider 
for the FLS field) and at the faintest 24/mi fluxes where the 
(random) differences rise to a ~ 0.5 mag at mi4 = 19 mag. 
We have performed aperture photometry in the Spitzer 3.6 
through 8 /mi IRAC images and 24/im MIPS images at the 
location of each of the SDSS/MPA galaxies, using a 7"radius 
aperture (12"at 24 /im). The IRAC fluxes are then aperture- 
corrected to 12.2", the radius used for IRAC standard star 
flux measurements. At 24 /im we attempt to calculate total 
fluxes, by applying an aperture correction to very large ra- 
dius (> 35") derived from the brightest 24/im sources. We 
have checked that the 7 "radius aperture (12"at 24/im) does 
not induce significant aperture effects by comparing to fluxes 
in larger radius apertures as a function of r band Petrosian 
radius; no significant trend is seen for the galaxies of our fi- 
nal sa mple. At 70/im our sources are unresolved, and fol- 
lowing [Gordoneral] d2007l) we apply an aperture correction 
of 2.07 to fluxes obtained with a 16"radius aperture with 18- 
39"sky annulus. Systematic errors in IR flux due to calibra- 
tion uncertainty, aperture corrections, and the resolved nature 
of many of the sources at wavelengths less than 70/im amount 
to - 30%. 

2.3. K-Corrections and Luminosities 

The resulting UV through 3.6/im magnitudes are K- 
corrected to z = 0.1, the median r edshift of the sample (e.g . 
V J g, etc.), using the method o flBlanton & Roweisl(|2007l) . 
The generalized UV luminosities are then calculated as L m , = 
AirDifxXf, where D^is the luminosity distance 14 and fx is 
the /T-corrected flux density per unit wavelength in the 01 f 
band, with effective wavelength A/ = 1390. At longer wave- 
lengths dust emission becomes more important than stellar 
emission, and we use a different method to ".g-correct" the 
data: we choose the best-fitting redshifted IDale et all (1200 ll) 
model IR spectral energy distribution (SED), on the basis of 
the observed 8 to 24/tm flux ratio, after correction for stellar 
emissi on at 8 and 24/im d eter mined from the 3. 6/im flux fol- 
lowing lHelou et al.ld2004 and lDale et al.l(l2005l) . This best-fit 
SED is then normalized using the measured 24/im flux, and 
the integrated far-IR (8-1000/i m) dust lum i nositie s (Ldust) are 
derived. Note that the different Dale et all (1200 ll) SEDs have 
LoAuxa/Ldust ratios that are different by a factor of up to five. 
We have checked that our results woul d not change s ignifi- 
cant ly if we use the model SEDs of iDevriendt et al.l (1 19991) 
(see iPapovich & Belll2002i for a detailed discussion of pre- 
dicting IR luminosities from Spitzer data). 

12 http://swire.ipac.caltech.edu/swire/astronomers.html for a description of 
the SWIRE image processing 

13 see http://ssc.spitzer.caltech.edu/postbcd/ for a description of post-BCD 
processing 

14 Throughout this study we assume a concordance cosmology: O,,, = 
0.3,n A =0.7,H = 70kms- 1 MpcT 1 . 



In Figure Q] we show the FIR flux ra tios for th e sam- 
ple galaxies, compared to the models of IDale et alJ (2001) 
and data from the Spitzer Infrared Nearby Galaxy Survey 
(SING S) sample of galaxies dKennicutt et alj2003l:rDale et all 
120051) . Note that the 24/im/70/im color does not vary as much 
as the 8/im/24/imcolor for these models and is not mono tonic. 
The 24/im/70/im ratio for the sample galaxies, and for the 
SINGS galaxies, sh o w a la rger range of L^(24/im)/L,,(70/im) 
than the IDale et al.l (2001) models. Many of the galaxies in 
our sample are bl uer in the 8/im/24/im color than any of the 
IDale et all d2001l) models at their redshift, likely due to in- 
complete treatment in these models of PAH emis sion, which 
is ass umed to be constant in spectral shape (see Smith et al. 
2007, for evidence of changes in the relative strength of PAH 
features). In general the scatter observed in this diagram sug- 
gests that mor e complicated and multi-parameter m odels may 
be necessary (Smith et al. 20071 iDraine & Lill2007l) . The use 
of such mod els with S DSS/MPA galaxies will be explored in a 
future work (I Johnson et al.l200"8T) . Some care must be taken in 
interpreting the total dust luminosities we derive, since these 
rely on models that do not appear to match existing data to 
high accuracy. 

2.4. Samples and Limits 

Our primary sample then consists of these galaxies with 
m 2 4 < 19.5(AB) and / and n < 25.0 (AB). It has been found 
that systematic errors in the n band photometry due to source 
blending can occur at n > 23 - only 45 of our galaxies are 
at these magnitudes. The SDSS spectroscopic selection, from 
which our sample is drawn, also requires r < 17.7 (AB). In 
the Lockman Hole, 80 of the 872 multiwavelength observed 
SDSS/MPA galaxies fail the UV cut, and 53 galaxies fail the 
24/im cut, of which 12 also fail the UV cut. We further require 
z > 0.02 (8 galaxies removed) and an r band Petrosian radius 
(R p (r)) of less than ll"(an additional 23 galaxies removed), 
to avoid large galaxies which are difficult to accurately pho- 
tometer. This leaves a total of 721 SDSS/MPA galaxies out 
of 872 observed at all wavelengths (/ through 24/im) in the 
Lockman Hole (82%). The galaxies that are missed due to 
the UV and IR flux cuts are primarily red sequence or early 
type galaxies (of which ~ 70% are detected), although sev- 
eral dwarf galaxies are also missed. The size cuts primarily 
remove blue dwarf galaxies - which tend to be at low-redshift 
in the SDSS sample - although a number of dwarfs remain. 
In the FLS these same limits yield 118 galaxies detected at 
all wavelengths out of of 186 observed (~ 63%). The lower 
detection rate in the FLS can be explained by the different 
UV exposure time distribution and the variation of the red- 
shift distribution between the FLS and Lockman Hole fields 
(which both have median redshift z = 0. 1). 

2.5. Classes of Galaxies 

We define several special classes of galaxies to be used 
throughout this study - these follow common definitions in 
the literature and are non-exclusive. 

• Dwa rfs, with logM* < 9, wh ere M* in solar units is 
from lKauffmann et al l (l2003bl) 

• Luminous infrared galaxies (LIRGs), with L<{ ust > 
10 U L Q 

• Ultraviolet luminous galaxies (UVLGs), with vL v = 
L uv > 1O 1O L . Note that this definition does not include 
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i/L M (24/Ltm)/i/L„(70/xm) 



FIG. 1. — FIR color-color diagram for the sample with 70/^m AB magnitude nvjou m < 15. Also sho wn are integrated values from Dale et al. (2005) SINGS 
galaxies (red triangles) and the models of Dale et al. (2001) (blue line, diamonds). The Dale et al. 1 2001) models are shown for redshifts of z — 0.0 (open circles), 
0. 1 (squares), 0.2 (diamonds), and 0.3 (stars). The median redshift of the sample is z = 0.1, Error bars show the photometric measurement errors. 



structural parameters (e.g. compactness) which ap- 
pear necessary in identifying galaxies that may serve as 
low redshift analogs f or Lyman break galaxies (LBG) 
dHeckman et alj|2005l) . 

We also consider the properties of galaxies that have been 
obser ved both as part of the SINGS program dKennicutt et al] 
120031) and by GALEX as par t of the Nearby Galaxy Survey 
(NGS: iGil de Paz et al.ll2007h . The in tegrated fluxes of these 
nearby, resolved galaxies are given by lDale et al. (2006]). The 
65 SINGS galaxies that we consi der have publicly avail- 
able drift-scan spectroscopy as well (Moustakas & KennicutH 
2006J, from which we have calculated a 4000A break mea- 
sure, D„(4000), for comparison with the SDSS D„(4000) val- 
ues (see 93. 11 1. 

3. THE GLOBAL PROPERTIES OF GALAXIES IN THE UV AND IR 

In lJohnson et al.1 (| 2006) we introduced an observational di- 
agnostic of galaxy properties that relates the SFH and atten- 
uation of a galaxy to its color. The SFH and attenuation are 
the primary drivers of galaxy color (although metallicity can 
play a significant role, especially in older galaxies). A combi- 
nation of IR and UV observations provides a measurement of 
obscured and unobscured SF, and is thus ideal for measuring 
SFR and attenuation. A detailed exploration of this relation 
is given in §|4] Here we investigate the location of nearby, 
resolved galaxies, special classes of galaxies, and the sample 
of all detected galaxies in the CMD and the relation between 
color, SFH, and attenuation. These additional diagnostic re- 
lations are the CMD and the relation between attenuation and 
SFR . The locations of the special classes of galaxies ( 92.5b 
both places them in the context of the larger, more general 
galaxy distribution and informs the interpretation of the rela- 
tions. The measures of SFH and dust attenuation that we use 
are the strength of the 4000 A break and the ratio of IR to UV 
luminosities, and are described here. 

3.1. D„(4000) 

As a measure of SFH we use t he definition of the 
4000A break strength given in Bal ogh et all (1 19991) 



(D„(4000)), since it is measurable at high signal to noise 
for the galaxies that we consider (i.e both red and blue se- 
quence galaxies), and is nearly insensitive to dust reddening 
(MacArthur 2005). The 4000A break is similar to a measure 
of the ratio of the SFR averaged over ~ 10 8 5 yr to the SFR 
averaged over > 10 9 yr but is subject to strong variations with 
metallicity, especially at larger break strengths. In $4] we ex- 
plore the behavior of D„(4000) for some simple galaxy star 
formation histories using stellar population synthesis models 
at two different metallicities. While there are some variations 
with metallicity, in general D„(4000) increases for smaller ra- 
tios of current to past averaged SFR (i.e "older" galaxies). 
One complication in the use of D„(4000) is that it is only mea- 
sured within the spectroscopic aperture, and so galaxies with 
moderate bulge/disk ratios may have D„(4000) overestimated 
- this is clearly t he case for a smal l numb er of galaxies in this 
sample (see also Kauff mann et al.l d2007l) ). 

3.2. IRX 

A common and robust attenuation indicator is the ratio of 
IR dust emission to UV stellar emission (the so-called infrared 
excess, IRX= log(Lj„ sf /L u ,,)). It is independent of and a qual- 
itatively different kind of measure than various optical mea- 
sures based on a color excess (e.g. the Balmer decrement). 
This indicator is predicated on the idea that the UV emission 
measures the amount of transmitted flux from young stars 
while the IR emission (here AA 8- 1000/im) measures the 
amount of UV flux from you ng stars absorbed by dust and 
re-radiated in the thermal IR dBuatHl992t ICalzetti et al]|1994t 
!KennicuttlfT998l: iGordon et alJl2000l)~ Under these simplistic 
assumptions the attenuation in the UV (in magnitudes) can 
be written A m , = 2.51og(10 w?z + 1). Light will be absorbed at 
wavelengths other than the UV, but the combination of dust 
attenuation curves that rise steeply towards the blue with blue 
intrinsic spectra has justified these assumptions when consid- 
ering intensely star-forming galaxies. For many of the galax- 
ies we consider this approximation may not be valid - specif- 
ically the fraction of the IR luminosity due to the heating of 
dust by UV light m ay vary between galaxies - but following 
iMeurer et al] d 19991) we define 
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A / ^=2.51og( ?? 10 IRX +l) 



(1) 



with i] = 1/1.68 for use when comparing the IRX we deter- 
mine to other measures of attenuation in magnitudes. The 
parameter 77 gives the fraction of the IR luminosity that 
comes from dust-absorbed UV light produced by young stars. 
Note that other conversion s between IRX and attenuation are 
available (e.g., iBelll 120021: iBuat et al.l l2005h but are gener- 
ally consistent with this formulation, which has the advan- 
tage of being physically motivated and simply understood. 
The value of ri that we choose is similar to that chosen by 



Buat et al. (2007) an d is consisten t with Meurer et al. ( 199' 
Belli d2003l).lHirashita etall (120031) . and llglesias-Ta ramo et al. 
(2006). In lJohnson et a l. (2007a) we attempt to estimate rj for 
each galaxy. 

In detail several modifications to the naive interpretation of 
IRX as a measure of attenuation are possible. First, there is 
the aforementioned effect of some fraction of the light being 
absorbed at wavelengths other than A ~ 1400A, which will 
vary depending on the SFH of the galaxy. Second, the possi- 
bility of different stellar populations within the galaxy being 
attenuated by differe nt amounts of dust must be considered 
(ICharlo t & Fall 2000). A third and related complication is the 
effect of the relative geometr y of the stars and the dust, on 
both l arge and s mall scales dWitt & Gordonl2000t IBuat & Xul 
119961) . although lGordon et al.1 d2000l) find that for young stel- 
lar populations IRX is a good measure of the UV attenuation 
for a range of dust geometries and extinction laws, but they 
do not consider the effects of inclinati on in realistic disks on 
this measure (e.g.. iPierini et al.ll2004l) . Finally, older, blue 
horizontal branch stars may contribute to the UV flux, and 
AGN may contribute to both the UV and IR emission. How- 
ever, despite these complications it appears that, at least for 
galaxies with D„(4000)< 1 .7, IRX is an accurate dust estima- 
tor that increases mo notonically with A f uv for a given SFH 
dJohnson et al.ll2007al) . 

It is important to clarify the meaning of IRX as used 
throughout this paper. For simplicity, and to conform with the 
convention in the literature, we will use IRX defined by the 
total IR and UV luminosities (IRX= Zog(L ( / UJ/ /L(/y), where 
L^asf and Lj/y are defined in §12. 3b . The issues mentioned 
above regarding the relation of IR X to the true FUV a tten- 
uation are considered separately in lJohnson et al.l d2007al) . 

3.3. Color-Magnitude Diagram 

The CMD for galaxies is a fundamental observational tool 
in the study of galaxy evolution. Neglecting the effect of 
dust (and to a lesser extent metallicity), the color coarsely 
measures the SFH of a galaxy, while the absolute magni- 
tude in one of the redder bands measures the total stellar 
mass (again subject to the effects of dust attenuation, and 
and to a lesser extent the SFH). The CMD thus relates the 
current star-f ormation to t he pas t-inte grated star formation 
of galaxies. iBaldry et al.l d2004 and iBlanton et all (120031) 
definitively identified the bimodality of galaxies in this di- 
agram at low redshift using the SDSS, with many galaxies 
havin g red or blue colors and few galaxies at intermediate 
color. iKauffmann et all (l2003bl) and lBrinchmann et all (120041) 
have constructed the SFH and stellar mass analog of the CMD 
for ~ 500,000 SDSS galaxies. Studies o f the evolution of 
the distribution o f galaxies in the CMD dBell et al.l l2004bt 
iFaber et aT1l2005l) have identified an increase in the mass den- 
sity of the red sequence from z ~ 0.7 to z ~ 0.1. 



FIG. 2. — The color magnitude diagram for the sample detected at all 
wavelengths. Top: The oi (n — 3.6/^m) color vs. M3 6. The color of each 
symbol encodes its IRX quintile, from lowest (purple) to highest (red). The 
median photometric errors are smaller than the symbol size. The inset shows 
the distribution of IRX and the bins used to color-code the symbols. Bottom: 
Same a s top, but with D„(4000) as the y-axis and M, from Kauffmann et al. 
1 2003b) as the jc-axis. Black diamonds show galaxies undetected in the UV 
or at 24pm, while squares mark galaxies with z < 0.02 or Petrosian radius 
greater than 1 l"in the r band. Error bars show the median error for the lowest 
quintile (purple) and highest quintile (red) of IRX, as well as the median error 
for the entire sample (black). 



In Figure [2] we show the CMD of the sample galaxies. 
In the top panel we use the n-3.6/im color - this is simi- 
lar to the traditional CMD s using u — r or g — r dBaldry et al.1 
l2004tlBlanton et aT]|2003l) . except that the use of the n magni- 
tude causes greater separation of the red and blue sequences 
dWyder et al.ll2007l) . The color of each symbol encodes its 
IRX quintile: red symbols are in the highest IRX quintile and 
tend to be located near the massive, high D„(4000) edge of 
the blue sequence or on the red sequence, while purple sym- 
bols are in the lowest IRX quintile and are found primarily 
in the low mass tail of the blue sequence (i.e., dwarfs). In 
the botom panel we use D„(4000) as the measure of SFH. 
This measure is more directly related to SFH than the broad- 
band color - as we will show below the color is additionally 
affect ed by dust attenuation. T he stellar mass estimates are 
from OCauffmannetiD d2003b) (ED- In this dia 

gram we 

also show the location of galaxies undetected in the UV or 
at 24/xm - these are approximately ~ 25% of the SDSS/MPA 
galaxies with D„(4000)> 1.7, as well as some dwarfs. Galax- 
ies that do not pass the size cuts (z > 0.02, R p (r) < 11") are 
also identified. 

We clearly see the red and blue sequences, and their sepa- 
ration, in both diagrams. The blue sequence is tilted such that 
it becomes redder with increasing stellar mass. The degree 
of tilt, the scatter around the main sequences, and the evo- 
lution of both can be used t o constrain the SFH of galaxies 
dSchweizer & Seitzerl 119921: IFaber et al.1 120051: iNoeske et al.l 
l2007t ISchiminovich et al.ll2007l) . There is also an increase 
in the average IRX with mass. Many of the galaxies in the 
'green valley' between the red and blue sequence appear to 
have large IRX. The accurate correction of observed colors 
for the presence of dust is thus a crucial step in determining 
the true tilt of the blue sequence for comparison to theories of 
galaxy formation. Attenuation is mass (and SFR) dependent, 
and so inaccurate attenuation corrections can lead to incor- 
rectly determined distributions of SFR with mass. 
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3.4. Dust-SFH-Color 

In Figure [3] we show the relation between IRX, D„(4000), 
and n-3.6/im color for the sample of galaxies presented in 
this work. The color of each symbol in dicates the IRX quin - 
tile. This diagram was introduced in iJohnson et al.l (2006). 
There is a clear relation between D„(4000) and color for 
the lowest quintiles of IRX (i.e., purple or blue symbols), 
though the relation becomes somewhat more scattered at high 
D„(4000). For larger values of IRX (i.e. orange or red points) 
the relationship between color and D„(4000) persists, but is 
shifted to redder colors because of the reddening effect of the 
dust. Figure [3] shows how the color is additionally affected by 
the attenuation as well as SFH - in this diagram the redden- 
ing vector is nearly horizontal - and how the conversion from 
color to D„(4000) in the CMD may be made, (i.e., to go from 
the top to the bottom panel in Figure |2j. 

The locations of LIRGs, UVLGs, and dwarfs are shown 
in Figure [3jz while the locations of nearby galaxies drawn 
from the SINGS and GALEX NGS are shown in Figure [2?. 
The LIRGs clearly occupy the the reddest region for a given 
D„(4000), and typically follow the sequence of galaxies with 
large IRX. The UVLGs are generally blue for their D„(4000), 
and have relatively low attenuation although they do not fol- 
low the sequence of galaxies with the lowest IRX. The nearby 
galaxies cover the parameter space probed by the current sam- 
ple reasonably well. 

In Figure [4] we show SDSS gri composite images of 
galaxies in the sample, selected randomly from bins of 
D„(4000) and 01 (n-3.6/xm) color. At high D„(4000) and rel- 
atively blue color we can see examples of galaxies with red 
bulges and blue, star forming discs, highli ghting the effect 
of ap erture on the D„(4000) measurement ( Kauffman n et al.1 
2007). At moderate D„(4000) and very red color there are 
many highly inclined, apparently red, discs, showing the ef- 
fect of attenuation on the color while D„(4000) is relatively 
unaffected. 

4. IMPLICATIONS OF THE RELATION BETWEEN ATTENUATION, 
SFH, AND COLOR 

Having identified the relation between dust attenuation, 
SFH, and color as a useful diagnostic of galaxy properties, 
here we seek to explore the nature of the relation in more 
detail by considering different projections and a wider range 
of colors sampling the entire UV-optical SED. We also com- 
pare the relation to models that predict, separately, the effects 
of dust attenuation and SFH on the color. In this way we 
can determine which UV-optical colors, in combination with 
D„(4000) and IRX, best constrain the models of attenuation 
and SFH. Finally, we parameterize the empirical relationship 
for these additional colors and examine the accuracy of the pa- 
rameterization, and the degree to which deviations are corre- 
lated with several additional physical properties of the galax- 
ies. 

4.1. Extension to Additional Colors 

In Figure we show D„(4000) versus color for a number 
of different broadband colors. The top panels show colors 
with a short wavelength separation, while the middle panels 
show colors with a longer wavelength separtion. As in Fig- 
ure [3] the symbols are color-coded by their IRX quintile. We 
see a clear separation of the effects of dust and SFH on the 
colors with widely separated wavelengths. This behavior is 
discernible in the colors with a shorter wavelength separation 
as well (e.g. g—r, u-g), but the separation is not as clear. For 



each panel of Figure|5]the effect of dust attenuation is to move 
galaxies nearly horizontally along the color axis, sinc e th e 
effect of attenuation on D„(4000) is minimal (but see 34.31 >. 
In effect, the 'reddening' vector is nearly horizontal. Thus a 
relation between D„(4000) and color that appears for galax- 
ies with low IRX (purple symbols) will simply be translated 
along the color axis for galaxies with large IRX (red sym- 
bols) if the relation between D„(4000) and unattenuated color 
is nearly constant among galaxies. The size of the translation 
is constrained by IRX. The clear separation of D„(4000) and 
dust effects also implies a nearly invariant relation between 
D„(4000) and unattenuated color, at least for regions of the 
diagram with low scatter (D„(4000)< 1.7). 

There are some changes in the form of the dust-SFH color 
relation depending on the color used. In particular, at red 
g — r color (g — r ~ 0.9) there is saturation in the sense that 
g—r color does not change significantly even as D„(4000) in- 
creases or for different values of IRX. This is in contrast to the 
n—r and other long-wavelength baseline colors where there 
is a relatively large variation in color for different values of 
D„(4000) and IRX. At low D„(4000) there is a difference in 
behavior of the n-r and n— 3 .6/im colors. The n— 3 .6/im color 
shows a larger range at low D„(4000) than the n-r color - this 
is because for a given D„(4000) and large dust attenuations the 
r band flux is affected more than the 3. 6/im flux, and so the 
n—r color does not increase as rapidly as the n — 3 .6/mi color 
for increasing attenuations. 

The trend of D„(4000) with f—n color for a given nar- 
row range of IRX is related to t he effect of SFH o n the UV 
spectral slope (3 as described by lKong et al.l (120041) . Here we 
see that there is very little, if any, correlation of f — n with 
D„(4000) for D„(4000)< 1.8 - bands of nearly constant IRX 
appear vertical in the D„(4000)-color plane, although there is 
a clear shift to redder colors for larger value s of IRX. The 
IRX- relation is considered in more detail in Johns onet al.l 
(I2007bl) . 

At D„(4000)> 1 .8 the relation between D„(4000) and color 
is different - there is a ridge of galaxies for which f—n de- 
creases for increasing D„(4000). IRX appears to be signif- 
icantly less correlated with f — n for these galaxi es tha n for 
those with D„(4000)< 1.8. As we see below ( 34.1.11 this 
ridge is close to the path followed by galaxies with a burst-like 
SFH and a small ratio of current to past averaged SFR, and 
traces the increasing importance of old blue stars in galaxies 
with very small ratios of young stars to old stars. 

4.1.1. SFH models 

To aid in the interpretation of the relation between 
D„(4000) and SFH - and between SFH and color - we have 
constructed several stellar p opulation synthesis ( SPS) m odels 
of galaxy spectra usi ng the iBruzual & Charlol (|2003) code. 
These models use the IChabrierl d2003l) IMF with the Padova 
1994 tracks, and have smooth, exponentially declining SFR^ 
e -'/ T sF w ith no gas recycling. Colors (in the blueshifted fil- 
ters) and D„(4000) are calculated directly from the result- 
ing spectra. No dust attenuation has been applied. We 
show four of these models in Figure [5] having metallicity 
Z/Zq = 0.4 or 2.5 and tsf = 0.1 or 1 Gyr. None of the mod- 
els are star-bursting in the sense that the ratio of the current 
to past-averaged SFR is always less than 1. A model with 
tsf = lOGyr follows tracks similar to the t$f = lGyr mod- 
els but never reaches D„(4000) > 1.4. In this diagram the 
models with tsf = 0.1 are similar to a simple stellar popula- 
tion (SSP) model. Models with different tsf generally occupy 
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FIG. 3. — Le/fc Relation between IRX, D„(4000), and color with location of LIRGs, UVLGs, and dwarf galaxies marked (see j|2.5l for definition of these 
classes). The color of each symbol indicates the IRX quintile, from lowest (purple) to highest (red) as in Figureff] Error bars in the lower right show the median 
eiTor for the lowest quintile (purple) and highest quintile (red) of IRX, as well as the median error for the entire sample (black). Right: Same as left, but with the 
location of galaxies from the SINGS and NGS surveys marked by grey diamonds. 



the same region of D„(4000)-color space at late times (t > 12 
Gyr), while models with different Z have significantly differ- 
ent D„(4000) at late times. It is important to keep in mind that 
the A"-correc tions we have appl i ed are derived from spectra 
based on the lBruzual & Charlotl ((2003) models, so some care 
must be taken in drawing conclusions based on small effects 
in color. 

The zero-dust SPS models nevertheless yield several inter- 
esting insights. First, the models with different t sf trace very 
different paths in the D„(4000) versus f—n diagram. For the 
tsf = 1 model the only change in f—n occurs at late times for 
Z/Zq = 0.4, even as D„(4000) steadily increases. This is con- 
sistent with the majority of galaxies at D„(4000)< 1.7 that, 
for a given small range in IRX, show little variation in f—n 
color even as D„(4000) varies. The tsf = 0.1 model exhibits 
strong variations in f—n color as a function of time. The 
models start at 01 (/-n) ~ 0.5 and become very red reaching 
a maximum of 01 (f—n) ^2 or ~ 2.5 for the high and low- 
metallicity models, respectively, at t ~ 1.5-2 Gyr. After this 
time the models become bluer as the emission from evolved 
blue stars becomes dominant, while D„(4000) continues to in- 
crease. Although several SDSS/MPA galaxies fall along the 
tsf =0.1 track at low D„(4000) and redder colors, these galax- 
ies also have significant IRX, such that they are incompatible 
with the zero-dust tsf = 0.1 model. The attenuation corrected 

(f — n) color of these galaxies would place them closer to 
the zero-dust tsf = 1 model. In contrast, many of the galax- 
ies with D„(4000)> 1.8 lie along a ridge that appears similar 
to the track of the fast declining [tsf = 0.1) model, and these 
galaxies appear to have little relation between their color and 
IRX. 

For the long wavelength separation colors ( 0A (n-r), 01 (n — 
z), and (n - 3.6/im)) the zero-dust tsf = 0.1 and tsf = 1 
models follow much more similar tracks in D„(4000)-color 
space. The galaxies in the lower IRX bins (purple and blue 
symbols) follow a similar track, although there is more scatter 
for D„(4000)> 1.7. The galaxies with the largest IRX (or- 
ange and red symbols) are redder than the unattenuated mod- 



els. Galaxies with significantly higher D„(4000) than any of 
the models for a given color are likely those for which the 
aperture effect causes an overestimate of global D w (4000), 
and h ave moderate bulge-to-disk ratios (e.g., Kauffma nn et alj 
120071) . The increased scatter at D„(4000)> 1.7 may be due 
in part to increased scatter in t he r elation between IRX and 
attenuation for such galaxies ( jj3.21 >. In addition, and impor- 
tantly, it is at D„(4000)> 1.7 where the models begin to di- 
verge, and our simple relation between D„(4000), attenuation, 
and color may be broken because of scatter in the relation be- 
tween D„(4000) and color. This scatter is in itself interesting, 
since it arises from different SFHs. It is only with robust, 
independent estimates of the attenuation that these small dif- 
ferences in the relation between color and D„(4000), corre- 
sponding to different forms of the SFH, can be determined 
from observations. The variation between the models of the 
relation between D„(4000) and unattenuated color for these 
longer wavelength separation colors suggests that these col- 
ors (as opposed to the optical colors, see below) give some 
clue to the details of the SFH. This makes the plots involving 
these colors, and the f—n color, the most useful for investi- 
gating the details of the SFH. 

For the colors 0A (u-g), 0l (g-r), and l (u-r) the be- 
havior of the models is different. These colors do not com- 
pletely straddle the redshifted 4000A break. For a given 
D„(4000) the different unattenuated models are nearly identi- 
cal in these optical colors, even at D„(4000)> 1.7. Because 
the expected, unattenuated relation between D„(4000) and 
color appears nearly independent of the details of the SFH 
for these colors, deviations from this relation are more reli- 
ably interpreted as due to dust attenuation. Thus, dust atten- 
uation would appear to be more easily inferred using the re- 
lation between D„(4000) and these colors than for the longer- 
wavelength baseline colors that sample the UV (e.g, 01 (n— z)). 
However, these optical colors also have larger errors relative 
to the length of the reddening vector than the long-baseline 
colors, which makes the inference of attenuation more diffi- 
cult. 
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FIG. 4. — SDSS gri color composite images of galaxies, organized in bins of D„ (4000) and 01 (n — 3.6/jm) color. For each bin of D,,(4000) or color (one eighth 
and one eleventh of the range, respectively) a galaxy is chosen at random for display. The top number in each image indicates the median IRX value for that bin, 
while the bottom number gives the number of galaxies in the bin. 



4.2. Binned by IRX 

In Figure [6] we show the relation between IRX and color 
for a number of different broadband colors. This is a different 
projection of the relation shown in Figures [3] and [5] and here 
the symbols are color-coded by D„(4000) quintiles. 

For the f — n color the relation bet ween IRX and color 
is analogous to the IRX-/3 relation dMeurer et al.1 119991: 
iKong et all |2004|) . where (3 is the power-law slope of the 
UV spectrum. The correlation between IRX and f — n has 
low scatter, and there is ver y little or no depe ndence on 
D„(4000) for D„(4000)< 1.7. IKong et al.1 (120041) have sug- 
gested that the additional effect of SFH on UV color may 
account for the scatter in the IRX-/3 relation. This would 
make the UV slope less reliable as a measure of attenuation 
in, e.g. high r e dshift galaxies, without knowledge of the SFH. 
iGordon etal.1 |2000) argue that the effects of geometry may 
make the UV slope a poor measure of attenuation. We find 
that, for a variety of galaxies with D„(4000)< 1 .7 observed in 
the local z < 0.3 universe, the f—n color is well correlated 
with IRX, although the slope is steep and therefore poorly 
constraining in the presence of significant errors in / - n or 
(3. However, for D„(4000)> 1.7 there is a large scatter in the 
relation, suggesting that the f—n color for these galaxies is 



not primarily driven by the attenuation ( 94.1. U or that IRX 
is not a good measur e of the attenuation for these galaxies 
(I Johnson et al.l l20~07al) . It is likely that both are true. The 
IRX-/3 r elation for this sample of galaxies is treated in more 
detail in lJohnson et al.l (l2007bl) . 

The g—r and u-g colors and near equivalents (such as U - 
V) are often used as i ndicators of red and blue sequence mem- 
bership in the CMD dBell et al.ll2004bl:lFaber et"al1l2005l) . We 
see here that they are somewhat affected by dust reddening 
dBell et al.ll2004al) . and that care must be taken in defining the 
red and blue sequence, especially if there is evolution in the 
distribution of attenuation of galaxies as a function of redshift. 

For the long wavelength separation colors (e.g. n - r or 
«-3.6//m) it is easy t o see the e f fect of SFH on the color 
that was predicted by Kong et al. (2004) for the UV color. 
For a given low D„(4000) (i.e. much recent SF, purple sym- 
bols) the relation between IRX and color is clear. This rela- 
tion is closely relat ed to the effective attenuation curve (e.g., 
iMeurer et al.fll999l) . As D„(4000) increases (i.e. for differ- 
ent colored symbols in Figure [6} the relation moves to redder 
color, but the slope and scatter do not change a great deal until 
the largest D„(4000). This shift to redder color at a given IRX 
is due to the redder intrinsic spectrum of a galaxy with larger 
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FIG. 5. — Relation between D„(4000) and color for different bins of IRX, for the colors / -n,u-g,g — r, ,n—r,n-3.6,f — r, u — r. The bins of IRX used to color 
code the symbols are shown at lower right. Error bars in the lower right of each panel show the median error for the low est quintile {purple) and highest quintile 
(red) of IRX, as well as the median error for the entire sample {black). Overplotted in black are four Bruzual & Chariot (2003) population synthesis models with 
Z/Zq = 0.4 (open symbols), 2.5 (filled symbols) and exponentially declining SFR with time constants t$f = 0.1 Gyr (circles), 1 Gyr (diamonds). The symbols 
are placed at intervals of 500 Myr, from 0.5 to 13 Gyr. 



D„(4000)(lower ratio of recent SFR to past averaged SFR). 
However, at the highest D„(4000) the relation between IRX 
and color becomes more scattered. 

4.2.1. Dust Models 

To more explicitly show the relation of Figure [6] to dust 
attenuation laws we have overplotted the relation between 
reddening and FUV attenuation given by the models of 
IWitt & Gordonl (l2000h . These models include several differ- 
ent local and global dust geometries, as well as both Milky 
Way (MW) and Small Magellanic Cloud (SMC) extinction 
laws. The several global dust geometries consist of equally 
mixed dust and stars in a sphere (DUSTY), a shell of dust 
surrounding a sphere of stars (SHELL), and a cloud of dust at 
the center of the more extended stellar distribution (CLOUDY, 
where the radius of the dust is 0.69 times the radius of the 
stellar sphere). The local distribution is described as either 
(h)omogenous or (c)lumpy. To place the model tracks in Fig- 
ure [6] we con vert the attenuation at ~ 1400 A to IRX using 
equation (1). iGordon et al.l (HOOO) show that this is a good 
approximation for all geometries with young stellar popula- 



tions, although as mentioned in $3.2l for the older galaxies in 
the sample this approximation may not be valid. We calculate 
the color excess by taking the difference of the attenuations at 
the effective wavelength of the filters. The zero point in color 
is determined by averaging the color of the five bluest galaxies 
in the lowest D„(4000) quintile. However, for "older" galaxies 
that will have a redder intrinsic, unattenuated color the model 
curves can be translated along the color axis to match the zero 
dust color for that galaxy - in other words, the horizontal ze- 
ropoint of these tracks is somewhat arbitrary, depending on 
the unattenuated spectrum of the galaxy. 

For clarity we only show the SHELL geometry in Figure 
[6] The CLOUDY model significantly underpredicts both IRX 
and the color excess for -ry < 10. While a few galaxies may 
be described by this model, the majority do no t appear to be in 
this co nfiguration, for the parameters used by I Witt & Gordonl 
(2000) (r dmt /r stars = 0.69). The same is true of the DUSTY 
model with a clumpy local geometry. The MW dust extinc- 
tion law is in significant conflict with the observed IRX versus 
f—n color relat ion. Because of the 2 175 A bump in the MW 
extinction law (IStecherl fl96l iDraine & Li 2007) the f-n 
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FIG. 6. — Relation between IRX and color for different bins of D„(4()00), for the colors f—n, u—g, g — r, n — r, n — z, n — 3.6/xm, f—r, and u — r. The bins 
of D„(40()0) used to color code the symbols are shown at lower right. Error bars in the lower right of each panel show the median error for the lowest quintile 
(purple) and highest quintile (red) of D n (4000), as well as the median error for the entire sample (black). The black symbols and dashed lines show the effect of 
dust attenuation from Witt & Gordon ( 2000), with a SHELL global geometry and for combinations of a MW (open symbols) or SMC (filled symbols) extinction 
law and a clumpy (squares) or homogeneous (circles) local dust distribution. Symbols run from Ty = 0.1 to 10 in steps of 0.5. See text for details. 



color becomes bluer with increasing Ty, a behaviour not ob- 
served in the data. As can be seen in the top left panel of 
Figure [6] and as has been noted by other authors, the rela- 
tion between UV slope and IRX is a very strong constraint 
on models of dust attenuation, particularly the extinction law. 
Note that some models that fit the data well for a given color 
do not fit the other colors. We find that the SMC extinction 
law with the SHELL geometry with clumpy or homogenous 
local geometries, or the DUSTY model with homogenous lo- 
cal geometry provide the best matches to the trends of IRX 
with color excess in the data, although the DUSTY model 
gives significantly lower IRX and color excess for a given Ty 
than the SHELL model. 

4.3. Binned by Color 

In Figure [7] we show the third projection of the dust-SFH- 
color relation, IRX versus D„(4000) for different ranges of 
n-3.6/im color. This effectively shows the region of IRX and 
D„(4000) parameter space allowed for a given color. The ex- 
tremely small amount of overlap of different color bins sug- 
gests that galaxies follow nearly parallel 'isocolor' lines in 



this diagram. Note that at high D„(4000) such isocolor lines 
do overlap, suggesting that an additional parameter may be 
necessary to explain the n-3.6^m color of these galaxies. 
This diagram also shows the correlation (or lack thereof) be- 
tween IRX and SFH in galaxies, and is comparable to similar 
diagra ms constructed f or resolved galaxies in the local uni- 
verse (Dal e et al.ll2006l) . However, one must again be cau- 
tious due to the effect of SFH on the interpretation of IRX as 
a measure of attenuation. 

The envelope of lowest D„(4000) galaxies moves to slightly 
higher D„(4000) at higher IRX, so that there are no galaxies 
with significant dust attenuation and extremely low D„(4000). 
It is possible that galaxies undergoing recent star-formation in 
the local universe rarely have much dust. However, it is also 
quantitatively consistent with the increase in D„(4000) ex- 
pected due to dust reddening, an indicati on that Dn(4000) i s 
not a perfectly 'clean' indicator of SFH dMac Arfhurl 2005 ) . 
Using an approximation for effective optical depth, Ty = 
[(5500/1390)" a7 /loge1-A /fi x where w e assu me r A oc A" 07 
we can compare a model of lMacAr thur (2005) to our sample 
in Figure [7] This is shown as a dashed line in Figure [7] This 
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FIG. 7.— Left: Relation between IRX and D„(4000) for different bins 
of (n — 3.6/xm) color. The dash ed line shows the e ffect of attenuation on 
D„(4000) for one of the models of lMacArthuJ 1200511 . Error bars in the top 
left corner show the median error for the lowest quintile of ■'(« - 3.6/^m) 
(purple) and highest quintile (red), as well as the median error for the en- 
tire sample (black). Right: Histogram of 01 (n — 3.6/^m) color, showing the 
quintiles used to color code the symbols at left. 

same trend can be seen in the data in Figure[5]where the enve- 
lope of low D„(4000) moves to larger values as the color (i.e. 
dust attenuation) increases. This shows the degree of redden- 
ing of D„(4000) for large values of attenuation, and gives the 
deviation of the reddening vector from an exactly horizontal 
line in this diagram. 

4.4. Parameterization of the Relation 

The relation between attenuation, color, and SFH described 
above is determined by the dust attenuation law and the cor- 
relation of D„(4000) with unreddened color. We leave a 
more rigorous derivat ion of a dust law to a separate work 
(I Johnson et al.1 l2007bl) . Here we give a simple parameteri- 
zation of the relation between IRX, color, and D„(4000) for 
a variety of colors. We fit using Amx as the measure of dust 
attenuation since this is linear in color, and should c lose ly ap- 
proximate the true Af uv for relatively blue galaxies Q3.2b . The 
relation between D„(4000) and color is allowed to be linear or 
quadratic, based on visual inspection of Figure[5] We are here 
assuming that the relation between D„(4000) and the unred- 
dened color is the same for all galaxies. Finally, after inspec- 
tion of the fit residuals, we allow a cross term whose origin is 
unclear, but may be related to the change in the relation be- 
tween Af uv and Amx as a function of SFH. Since the errors 
in Amx are larger than those in color or D„(4000), we treat 
Ajrx as the dependent variable. Thus, we fit 

A !RX = A+Bx + Cy+Dx 2 + Exy (2) 

where x =D„(4000)-1 .25 and y is the AB color, minus 2. The 
A term thus provides the typical attenuation for galaxies in 
the sample with D„(4000)= 1.25 and a color of ~ 2. The re- 
sults of these fits are given in Table [2] for several different 
colors. We have also conducted fits after restricting the sam- 
ple to D„(4000)< 1.6 and 01 (n-r) < 4. This allows us to 

ignore the effect of 'older' galaxies for which Amx is likely 
not a good indicator of attenuation, and for which dust red- 
dening may play a smaller role in determining the color, than, 
e.g. metallicity, AGN, or blue evolved stars. The inclusion of 
these high D„(4000) galaxies may bias the fits, and when con- 
sidering only blue sequence galaxies we encourage the reader 
to use the fits restricted to D„(4000)< 1.6. 

The fits including a cross term (E) but without a quadratic 
term in D„(4000) (D= 0) are shown in Figure [8] along with 
the sample galaxies. This figure shows fits for both the en- 



tire sample and for the D„(4000)< 1 .6 subsample, in all three 
projections of the relation, and for the four colors given in Ta- 
ble |2 Because the displayed fits do not include D term, the 
curvature seen in all colors is due to the cross-term, although 
for '("-3.6/mi) the fit to galaxies with D„(4000)< 1.6 ap- 
pears very linear. The cross-term also causes the lines for 
different Amx values to not be parallel. Despite the compa- 
rable dispersion in the residuals, the fits using other combi- 
nations of terms are not as well behaved. Residuals of the 
fits to the entire sample of galaxies are shown as a function 
of D„(4000) and color in Figure [9] It is possible that addi- 
tional parameters play a role in the relation. For example, the 
metallicity of galaxies may have a significant additional ef- 
fect on the color beyond the effects of SFH and attenuation. 
As a preliminary investigation of this possibility we show the 
residuals from the fit using ()1 (n-3.6^m) plotted against other 
galaxy parameters in Figure[l0] Th ese parameters are th e gas- 
phase metallicity as determined bylTremonti et alJ (120041) . the 
stellar mass from Kauffmann et al. (2003b), the SFR as deter- 
mined by Brinchman n et alJ (12004) (where we use the median 
of the probability distribution function for the global SFR), 
and A/flx-The strongest trend is with Amx, the very quantity 
being fit. This is a result of the dispersion in the relation being 
a large fraction of the range, and precludes detailed analysis 
of the residuals. 

4.5. Applications and Implications for Future Observations 
and Models 

The A/Rx-color-D„(4000) relation may be used to determine 
the attenuation given D„(4000) and a color. Such derived at- 
tenuations may then be used to correct SFR indicators at other 
wavelengths. While D„(4000) is already a coarse measure of 
the star-formation activity of galaxies, it is sensitive to older 
stars than the / band or Ha luminosity. The accuracy of 
the attenuation correction derived from the fits and applied to 
UV or Ha luminosities then depends on whether the most re- 
cently formed stars are affected by the same dust as those stars 
contributing to the long wavelength separation color (as well 
as the accuracy of the fit for the chosen color, and whether 
the form of the relation between D„(40 00) and color varies 
with changes in the details of the SFH - ). iTreyer et al] (2007) 
have compared the UV fluxes corrected for attenuation using 
these fits to the Ha derived SFR. They find good agreement, 
although there is a trend of the differences with the SFR. A di- 
rect comparison of the UV+IR derived SFR to the emission- 
line de rived SFR is given in a separate work dJohnson et alJ 
l2007ah . 

Here we compare the accuracy of this method for attenu- 
ation correction to methods based on emission lines or the 
UV color. The standard deviation of the residuals from the 
fits, in terms of both Amx and IRX, are given in Table [2] 
For the 01 (n - 3.6/im) color these are 0.5 in Amx and 0.3 
in IRX. Attenuation may also be estimated from the Balmer 

decrement A Ba i mer = 2.5 log where L Ha and h H s 

are the Ha and hb line luminosities, repectively dKennicuttl 

1998; Johnson et al. 2007a). The ratio AmxlA Ba imer is ~ 3-5 
dJohnson et alJ l2007ah . Thus, to obtain a similar error in 
Amx from the Balmer decrement requires S/N > 10 in both 
the Ha and H/3 emission lines (or an equivalent combina- 
tion of S/N in the two lines). Another common method of at- 
tenuation correctio n is provided by the Af uv -f3 relation (e.g., 
iMeurer et alJI 19991) . where (3 is the UV spectral slope (simi- 
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FIG. 8. — Fits of Table |2]overplotted on the 3 projections of the relation between IRX, D„(4000), and color. From left to right the colors used are 01 (n — r), 
01 (n— z), 01 (n-3.6fim), o i (g—r). The solid lines show the fits determined from galaxies with D„(4000)< 1.6 and 1 (n — r) < 4, while the dashed lines show the 
fits determined from the entire sample. 



lar to the f—n color). ISeibert et alJ (120051) quote a scatter of 
a(Af uv ) = 0.89 in this relation when a range of galaxy types is 
considered. One must also consider the large slope in the rela- 
tion between (3 and A/,,,,. An uncertainty of 0.1 in/— n color 
(due to photometric errors or uncertainty in the /^-correction) 
leads to an uncertainty of ~ 0.8 in Af uv , even in the absence 
of scatter in the Af uv -j3 relation. Such uncertainties in color 

(and D„(4000)) are much less important for the A/sx-color- 
D„(4000) relation, especially when UV-optical colors are used 
that have large wavelength separation and relatively small er- 
rors. 

Besides correcting SFR measures at various wavelengths 
for dust attenuation, it may be interesting to consider the dis- 
tribution of dust attenuations themselves - derived from the 
A/Rx-color-D„(4000) relation - as a probe of the star forma- 
tion in galaxies. The attenuation may be expect ed to trace 
the product of metall icity and gas surface density) Belll l2003L 
iJohnson et alJl2007bh . whic h are both key di agnostics of the 
evolution of galaxies (e.g., iMartin et al.l 120061) . The fit pre- 
sented in 34.41 can be used to estimate the dust attenuation 
for large samples of galaxies, at a range of redshifts, where 
D„(4000) and a broadband color are known - high S/N emis- 
sion lines are not necessary. 

Several future observatories will add to the IR view of 
the univ erse provided by S pitzer, including AKARI (formerly 
Astro-F, iMurakarnl 1 1 9981) and the proposed Wide- Field In- 
frared Survey Explorer (WISE; iMainzer et al.ll200"5l) . In par- 



ticular, the flux limit of WISE at 23/im (~ 2.6 mJy, 5a) will 
enable the detection of ~ 20% of SDSS/MPA galaxies at this 
wavelength (i.e., 10 5 galaxies), although these will be pri- 
marily dusty star-forming galaxies. The parameterization of 
the relation given above will provide a context for this much 
larger sample of galaxies that will be strongly affected by se- 
lection effects in the IR. Combined with GALEX imaging over 
the SDSS footprint, the large number of galaxies can be used 
to investigate deviations from the relation as a function of, 
e.g. AGN activity or morphology - understanding such de- 
viations will lead to improved modelling of attenuation and 
knowledge of the SFH. Such a large number of galaxies will 
also provide an unbiased view of the properties of rare types 
of galaxies (e.g., ULIRGS), which may or may not follow the 
A/Rx-color-D„(4000) relation derived here for more common 
galaxies. LIRGs at z = 0.5 will be detectable by WISE over the 
entire sky at 23fim, while the 12/im band of WISE will be able 
to probe the restframe 8/im PAH emission of these galaxies. 
Furthermore, nearly all SDSS/MPA galaxies will be detected 
in the shortest wavelength band (~ 3.3/im). Such NIR pho- 
tometry - when combined with resframe near-UV photome- 
try - leads to the clearest relation between broadband color, 
D„(4000), and IRX, making accurate determinations of IRX 
possible for ~ 10 5 galaxies. 

Finally, this parameterized relation between D„(4000), 
IRX, and color will have to be reproduced by the next gen- 
eration of galaxy models that seek to self-consistently treat 
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FIG. 9. — Residuals from the fit (Table [3) for each color, plotted vs 
D„(4000) and color. The residuals shown are from the fit to the entire sample 
including the A, B, C, and E terms in Table [2] Dashed lines show the stan- 
dard deviation of the residuals. The colors of the symbols in the left panels 
give the quintile of the observed color, while the color of symbols in the right 
panels gives the D„(4000) quintile (see Figure[6] bottom right panel). 

dust absorption and emission in the context o f SPS modeling 
dJonssonl2006l:lLi et alj2007l:ISilva et aljl998l) . The emprical 
relation given above provides an important diagnostic for such 
models since it links, in a relatively straightforward way, the 
stellar populations, the absorption of the light of these popu- 
lations by dust, and the subsequent re-emission of that light 
by dust at IR wavelengths. 

5. COLOR-COLOR RELATIONS AT HIGH REDSHIFT 

At high redshift (z ~ 2) several authors have attempted to 
distinguish between the effects of attenuation and SFH on the 
colors of g alaxies using Spitzer IRAC imaging of the rest- 
frame NIR (IWuyts et alJl2007t iKriek et alJl2006l) These use 
a methodology similar to that of iKauffmann et all (l2003bl) by 
comparing the locations of galaxies in color-color space to a 
locus o f SPS models, and attributing differences to dust red- 
dening. IWuvts et~aT] (12007b have used the relatively reddening 
insensitive photometric break strength at 4000A (sampling 
both D„(4000) and the Balmer break) explicitly as one of the 
'colors', with the benefit that the direction of the reddening 
vector is then relatively well known. The determination of 
IRX, as in the present work, gives empirical constraints on the 
length of the reddening vector, lending support to this method 
for SFH determination. 

In Figure QT| we show a color-color diagram for the sam- 
ple galaxies that is quite similar to that constructed by 
Wuyts ^tail (120071) for a number of galaxies at z = 1-3, both 
in the colors plotted, and in the location of galaxies in the 
diagram. Galaxies tend to move toward redder colors in 
01 (g- 3.6/im) as IRX increases, and are generally at high 
D„(4000) as ol (u-g) increases. As we have shown previously 
(" 94.1.1b there is also some reddening of the u — g color with 



IRX. For a large portion of the color space (i.e. 01 (w — g) < 
1.5), it is not possible use the color-color diagram alone to 
distinguish between different model tracks and the amount of 
dust attenuation. With independent constraints on dust atten- 
uation, and assuming that the population synthesis modeling 
yields the correct colors, it is possible to rule out high metal- 
licity (Z/Z = 2.5) exponentially declining models for a large 
proportion of the sample, although low metallicity exponen- 
tially declining models are consistent with the data. 

6. CONCLUSIONS 

1. We have presented a low redshift (z= 0.1) context for 
the IR, UV, and optically derived dust attenuation, SF, 
and stellar mass properties of ~ 1000 homogeneously 
observed galaxies. The relation between SFH, attenua- 
tion, and color is a useful diagnostic of galaxy proper- 
ties. We identified where special classes of galaxies lie 
in the "color-color" diagram of D„(4000) versus broad- 
band color. We find that LIRGs occupy the same region 
of D„(4000), IRX, and color parameter space as dusty 
star-forming galaxies, of which many are likely to be 
highly inclined. We find that UVLGs occupy the same 
parameter space as a larger population of relatively at- 
tenuation free, star forming galaxies. However, UVLGs 
are not the least dusty galaxies for a given D„(4000). 
The distribution of galaxies in this parameter space, 
corrected for selection effects, can be compared to high 
redshift samples and semi-analytic models to determine 
the evolution of dust attenuation, and to insure that at- 
tenuation is being treated correctly in the models. 

2. This relation between D„(4000), IRX, and color sug- 
gests that the colors of galaxies can be explained pri- 
marily as the contributions from the star formation 
history (we have here assumed expone ntially declin- 
ing m odels) and the dust attenuation dJohnson et al.l 
2006). A comparison with stellar population synthe- 
sis and dust attenuation models confirms that varia- 
tion in D„(4000) and IRX causes changes in galaxy 
color quantitatively similar to those seen in the data. 
Detailed comparison of the exact form of the relation 
of D„(4000) and IRX to the color of galaxies has the 
power to constrain these models. The UV color is es- 
pecially sensistive to the for m of the extinction law - 
a Milky Way extinction law dCardelh et alj|1989t) is not 
consistent with the change in UV colors of galaxies as a 
function of I RX, for a range of rel ative star/dust geome- 
tries (but see lPanuzzo et al.|[2007l for a possible resolu- 
tion involving the variation of attenuation with stellar 
age in galaxies). 

For some colors (e.g. 01 (u-r)) the relation between 
D„(4000) and color appears to be independent of the 
specific SFH model assumed. Deviations of the ob- 
served galaxy colors from this relation may then be ex- 
plained as reddening due to dust attenuation. However, 
for other UV-optical-NIR colors (e.g. oi ( w_ 3.6/jm)) 
the relation between D„(4000) and color depends on 
the specific form of the SFH, especially at moderate 
D„(4000) and color. In this case, independent con- 
straints on the degree of reddening provided by IRX 
will allow different SFH to be distinguished. Similarly, 
different SFH models follow extremely different tracks 
in the D„(4000) versus UV color plane, emphasizing 
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FIG. 10. — Residuals from the fit to the full sample, including the A, B, C, and E terms in Table|2] plotted against other galaxy parameters. The colors used 
in the fit are (top to bottom), 01 (« — r), 01 (n — z), 01 (n— 3.6/um), and (g—r). The residuals are plotted against (left to right), the gas-phase metallicity (where 
available), the stellar mass, the optically determined SFR (from Brinchmann et al. 2004), and Airx. 

parametrization may be used to infer the distribution of 
IRX in galaxies, and the evolution thereof from high 
to low redshift. It may also be used to determine the 
dust attenuations for application to tracers of more re- 
cent star formation. The best accuracy in inferred at- 
tenuations is achieved with the 01 (n-3.6/im) color - 
AKARI and the proposed WISE observatory, combined 
with GALEX will provide large samples of galaxies for 
which this color is available. This relation links dust 
absorption and emission and the stellar populations of 
galaxies. Models that seek to explain the optical and in- 
frared emission of galaxies self-consistently will need 
to reproduce this relation for low-redshift galaxies. 
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FIG. 11. — Color-colo r diagrams as are available at high redshift using 
IRAC <Wuvts et al] f20071. On the left the symbols are color-coded by the 
D„(4000) quintile, on the right they are color coded by IRX quintile . See 
Figures [6] and \5\ for the bins used for color coding, respectively. The model 
symbols (black) are the same as in Figure \5\ Colored error bars at the top left 
corner of each panel show the median errors for the corresponding parameter 
bin, while the black error bars show the median for the entire sample. 

the power of the UV color as a diagnostic of stellar pop- 
ulations. 

3. We find through comparison with zero-dust stellar pop- 
ulation synthesis models that a burst-like SFH leads 
to a relation between D„(4000) and UV color at late 
times that is similar to that seen for our sample galax- 
ies at D„(4000)> 1.7. The colors of these galaxies 
do not appear to be highly correlated with IRX, sug- 
gesting that dust attenuation is not the dominant driver 
of UV color for these galaxies. This effect is likely 
caused by evolved blue stars. There is, however a sig- 
nificant offset of the locu s of sample galaxies from the 
iBruzual & Charloil (12003) models of old stellar popula- 
tions. 

4. We have parameterized the relation between D„(4000), 
color, and IRX for a variety of colors. This allows 
the determination of one of these quantities provided 
the other two are known. Assuming a universal rela- 
tion between D„(4000) and unattenuatted color (as pro- 
vided by most exponentially declining SF models), this 
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TABLE 2 

Empirical Fits for Different Broadband Colors 



All Galaxies D„(4000)< 1.6 and " '(n-r) < 4 (401 galaxies) 

A B C D E <r(A, RX ) cr(IRX) I A B C D E a(A IRX ) ct(IRX) 















>.=°'(n 


-r)-2 














1.61 


-2.96 


0.77 


— 


— 


0.76 


0.40 


1.21 


-3.16 


1.38 


— 


— 


0.60 


0.33 


1.48 


-1.12 


0.76 


-2.50 




0.74 


0.38 


1.24 


-2.12 


1.36 


-4.65 




0.57 


0.32 


1.25 


-1.33 


1.19 




-1.02 


0.70 


0.35 


1.20 


-2.51 


1.48 




-1.05 


0.59 


0.32 


1.20 


-3.35 


1.57 


4.69 


-1.91 


0.68 


0.35 


1.22 


-2.23 


1.44 


-2.00 


-0.79 


0.59 


0.32 














y= 0l (n 


-z)-2 














1.23 


-3.27 


0.78 






0.73 


0.38 


0.74 


-3.26 


1.18 






0.58 


0.32 


1.13 


-1.66 


0.76 


-2.14 




0.72 


0.37 


0.77 


-2.28 


1.16 


-4.34 




0.57 


0.32 


0.81 


-1.32 


1.07 




-0.81 


0.68 


0.35 


0.71 


-2.49 


1.24 




-0.67 


0.57 


0.32 


0.65 


-2.79 


1.36 


4.22 


-1.53 


0.66 


0.35 


0.75 


-2.20 


1.20 


-2.90 


-0.35 


0.57 


0.32 














y= 01 (n- 


3.6jum)- 


2 












1.25 


-3.49 


0.84 






0.63 


0.35 


0.98 


-2.92 


1.03 






0.46 


0.28 


1.17 


-2.24 


0.83 


-1.64 




0.62 


0.32 


1.00 


-2.28 


1.02 


-2.87 




0.46 


0.27 


0.98 


-1.89 


1.02 




-0.63 


0.59 


0.31 


0.98 


-2.65 


1.05 




-0.25 


0.46 


0.27 


0.94 


-2.80 


1.16 


2.51 


-1.02 


0.58 


0.31 


1.00 


-2.27 


1.02 


-2.81 


-0.02 


0.46 


0.27 














y= 0l ( 


g-r)-2 














10.23 


-3.17 


5.75 






0.74 


0.38 


10.60 


-3.71 


5.99 






0.67 


0.36 


10.47 


-3.63 


5.90 


0.51 




0.74 


0.38 


10.46 


-2.85 


5.89 


-3.67 




0.67 


0.36 


10.35 


-4.34 


5.85 




-1.05 


0.74 


0.38 


10.84 


-6.64 


6.15 




-2.10 


0.67 


0.36 


12.05 


-10.59 


6.96 


2.76 


-4.44 


0.74 


0.38 


10.38 


-1.98 


5.84 


-4.18 


0.54 


0.67 


0.36 



